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We have investigated the dynamics of underdamped Josephson junctions. In addition to the 
usual crossover between macroscopic quantum tunnelling and thermally activated (TA) behaviour 
we observe in our samples with relatively small Josephson coupling Ej, for the first time, the 
transition from TA behaviour to underdamped phase diffusion. Above the crossover temperature 
the threshold for switching into the finite voltage state becomes extremely sharp. We propose a 
(T, Ej) phase-diagram with various regimes and show that for a proper description of it dissipation 
and level quantization in a metastable well are crucial. 
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A hysteretic Josephson junction (JJ) switching from 
its metastable zero-voltage state into a stable state with 
a voltage of the order of twice the superconducting gap 
A has recently been used as a read-out device for super- 
conducting quantum bit systems in many experiments 
Switching measurements also enable one to perform 
conventional large bandwidth current measurements, and 
recently there have been proposals to use hysteretic JJs 
as ammeters for studying phenomena like non-Gaussian 
noise 2J. Often it may be advantageous to reduce the 
critical current Ic of the detecting junction in order to 
increase the measurement sensivity. Yet the physics gov- 
erning switching phenomena of small Ic junctions ulti- 
mately differs from those with larger 7^, and this sets a 
limit on how far one can reduce L- still maintaining the 
useful features of the detector [3, ■ 

The behavior of a JJ is determined by various param- 
eters. Within the RCSJ (Resistively and Capacitively 
Shunted Junction) model these are the Josephson energy 
Ej = hlc/2e, the charging energy Ec = e^/2Cj (Cj 
is the junction capacitance), and a shunt resistance R 
responsible for dissipation. The dynamics of the junc- 
tion are then characterized by the plasma frequency 
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y/8EjEc/h and the quality factor Q 



As we will detail below, the behavior of the junction at 
a given temperature T can be classified according to the 
phase diagram of Fig. The overdamped case Q < 1 
was studied in detail by Vion et al. 0], who uncovered 
the existence of a phase diffusion regime at finite T with 
the appearance of a small voltage, prior to switching 
to a voltage of the order of 2 A. As far as the under- 
damped case Q > 1 is concerned, most experiments were 
done on large hysteretic junctions with relatively high Ic- 
Depending on temperature such junctions escape out of 



the metastable zero-voltage state either via macroscopic 
quantum tunnelling (MQT) or thermal activation (TA) 
processes, thereby switching directly to the finite- voltage 
state. In this Letter we focus on the regime Q > 1 with 
relatively small Ic- We show, theoretically and exper- 
imentally, that a regime exists where escape does not 
lead to a finite voltage state but rather to underdamped 
phase diffusion (UPD, the shaded region in Fig.QJ. 

According to the RCSJ model, the dynamics of a JJ 
biased with a current / is that of a phase particle (phase 
ip) with mass m — fi^ /8Ec, moving in a tilted cosine po- 
tential U{(p) = —Ej {cos(p+I / Ic (f) under a viscous force 
(^)^ 1/R The bias current renormalizes the plasma 

frequency, such that ujp = \f^^^ = j-^/SEjE^ q^^"^ 

with qo = ■y/2(l — I /Ic)- The cosine potential has wells 
where the phase particle can be localized; the phase 
then has constant average value, and the average volt- 
age across the junction is zero. For non-zero / the 
quantum levels in the potential well are metastable and 
the particle can escape from a given well either via TA 
over, or MQT through the barrier. For large junctions 

Tta = 2^ ^ ''B^ for the TA escape rate and Tmqt = 

, ,, , 36 AiJ 1—1 

\2^/&^'^^MJ/nwpe ^ ""p for the MQT rate 0. In the 
cubic approximation the barrier height is AC/ = | Ej q^. 
Below the cross-over temperature Tq — tiujp/2'KkB the 
dominant escape mechanism is MQT. The total escape 
rate is r(/) ~ Tmqt {I) + ^ta{I) and the total escape 
probability in the time interval < < < r can be written 
as P{I) = l~e^ ^0 r(/)rft^ If dissipation is weak, upon es- 
cape from the well the particle moves down the potential 
and phase is running freely, hence the voltage reaches a 
finite value (about 2A^57g« 360 /uV). However, if dissi- 
pation is strong enough, there is a finite probability that. 
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FIG. 1: (a) The various operation regimes of a Josephson 
junction with low Ej. For details, see text, (b) The cosine 
potential with dynamics inside the upper well and schematic 
dynamics after leaving the upper well. Inset: equivalent cir- 
cuit of the junction with frequency dependent dissipation. 



upon escape from the well, the particle is relocalized in 
the next well instead of running down the potential: the 
phase then has 27r-slips and diffusively moves from one 
metastable well to another. In this UPD regime the av- 
erage voltage across the junction is much smaller than 
2 All. 

A misconception still persists that phase diffusion is 
absent in a hysteretic junction, although it was pointed 
out over a decade ago that it can occur due to the depen- 
dence of dissipation on frequency lo 6] . Our experiment 
corresponds to the simplified equivalent circuit with fre- 
quency dependent dissipation as presented in the inset 
of Fig. nib). After switching to the running state, the 
dominant part of dissipation comes from small lo, gov- 
erned by R{lu ^ 0), typically given by the large junction 
subgap resistance, of the order of 1 M^l. In the phase 
diffusion regime the phase mainly oscillates in a well at 
the plasma frequency and thus the dissipation is charac- 
terized by Riujp), which is much smaller, typically of the 
order of vacuum impedance Zq w 377 fi, since Cg acts as 
a short. In the following we will consider junctions that 
are underdamped even at LOp, in contrast to 3. 

The dissipated energy between neighboring potential 



maxima can be approximated by Ed « 8Ej/Q and if 
the phase particle has energy less than above the 
next barrier top, it simply diffuses to the next well. The 
maximum possible dissipated power due to phase diffu- 
sion can be written as -^^^Eu, where V is the average 
voltage across the junction. By equating this with the 
applied bias power ImV, we find the maximum possible 
phase diffusion current 



(1) 



which is identical in form to the well known retrapping 
current formula, but now the value of Q is that at plasma 
frequency ujp. For I < Im, there is non-zero probability 
that phase relocalizes after escape. With decreasing /c, 
dissipation starts to play a more significant role. First, 
ujp and hence the quality factor Q are decreasing with 
decreasing Ej. On the other hand the escape rate is sig- 
nificant in the range of currents where the barrier height 
is comparable to the thermal energy. For large junctions, 
E,j ^ ksT, this is the case when / is just slightly below 
Jc, but for small junctions this occurs even without tilt- 
ing the cosine potential {I/Ic=0)- At small values of Ic 
the dissipation after escape is thus also important, be- 
cause the successive barrier tops are close in energy. The 
grey area in Fig. ^a) presents the UPD regime, where 
escape does not necessarily lead to the transition into a 
running state. The minimum Ef, where such a transi- 
tion after escape is still certain at temperature > Tq 
with current pulses of length r. 



E' 
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determines the line separating TA and UPD regions in 
Fig. n^a). For T < Tq, the escape rate equals the TA rate 
at To, thus Ef is independent of T, and given by Eq. Q 
with To replacing Tu- The UPD region in Fig. ^a) cor- 
responds to Rs= 500 n, Cj= 100 fF and T = 100 us. 

We present experimental data of two samples, a dc- 
SQUID and a single JJ. The dc-SQUID consists of two 
wide superconducting planes connected by two short su- 
perconducting lines with tunnel junctions in the middle 
forming the dc-SQUID loop of area 20 x 39 (/xm)^ (see 
the inset in Fig. [2J|. The loop inductance was measured 
to be around 100 pH, small as compared to the calculated 
Josephson inductance {L = = 400 pH per junc- 

tion). The dc-SQUID thus behaves almost like a single 
JJ, whose Ic can be tuned. The loop inductances were es- 
timated from the measured resonant voltage determined 
by Cj and loop inductance 0] . The other measured sam- 
ple was a single junction between long inductive bias- 
ing lines. The normal state resistances of the dc-SQUID 
and the single J J were 1.3 kfi and 0.41 kO yielding for Ic 
199 nA and 630 nA, respectively. The geometrical capac- 
itances of the samples were 100 fF and 130 fF. The mea- 
sured samples were fabricated using standard electron 
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beam lithography and aluminum metallization in a UHV 
evaporator. Both measured samples had strongly hys- 
teretic /1^-characteristics with retrapping currents well 
below 1 nA. 

The experimental setup is presented in the inset in Fig. 
|5| Switching probabilities have been measured by apply- 
ing a set of trapezoidal current pulses through the sample 
and by measuring the number of resulting voltage pulses. 
At the sample stage we used low pass i?C-filters (surface 
mount components near the sample). In the measure- 
ments on a single junction we used surface mount capac- 
itors (Cs = 680 pF), but in the dc-SQUID measurements 
we had 7r-filters in series with resistors , with ~ 5 nF 
capacitance to ground. The resistors were Rs — 500 ft 
and 681 fl in the measurement on a dc-SQUID and on a 
single junction, respectively. Bonding wires with induc- 
tance of order nH connect the sample to filters. 
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FIG. 2: Cumulative histograms of the dc-SQUID at differ- 
ent temperatures and at different magnetic fields. Curves 
are shifted for clarity and the vertical spacing between ticks 
corresponds to unity escape probability. Solid lines are from 
simulations described in the text; dotted blue lines show the 
results of basic MQT and TA models. Inset: scanning electron 
micrograph of the measured dc-SQUID and the experimental 
circuit. 



In Fig.|21we present the measured cumulative switching 
histograms (open circles) of the dc-SQUID at different 
fluxes $ and temperatures with r — 200 /is. At the low- 
est temperatures, the histograms can be well fitted by the 
MQT model, giving 4 = 200 nA, 128 nA and 55 nA, for 
<i>/<i>o = 0, 0.28 and 0.41 respectively. For $ = wc also 
present the simulated histograms based on the basic TA 
and MQT models without dissipation (dotted blue hues). 
In Fig. El we show the measured histogram position I^q% 
[Piho%) = 0.5] and the width A/ (= - /io%) for 
both samples. The dc-SQUID measurements were done 



both at negative and positive values of flux in order to 
ensure that the external flux had not changed. The posi- 
tion in Fig. Ofa) is normalized to the corresponding value 
of Ic at zero temperature. We also present the width and 
position results of the basic TA and MQT model simula- 
tions. At low T all the measured data are consistent with 
MQT results. On increasing T the parameters are con- 
stant up to the estimated cross-over temperature Tq. For 
T > Tq the width is increasing and the position is moving 
down as TA model predicts. The qualitative agreement 
is good for most of the results up to the temperature 
Td. At Tu, A/ starts to decrease abruptly. Moreover, 
the position 1^0% saturates at the same value ~ 0.35/c. 
For the single JJ, the saturation occurs at ~ 0.3/c. The 
dc-SQUID data measured at $ = ±0.41$o are not fol- 
lowing the standard theory even at low temperatures, 
since Tjj < 30 mK in this case. If we assume a realistic 
shunt impedance R{ujp) ~ 500 Q (the value of the sur- 
face mount resistors) , we obtain Q « 4 (4 = 200 nA 
and Cj = 100 fF), which yields /„,. ~ 0.35/c through 
like in the experiment. In the diagram of Fig. ^ we 
also present Ic of the dc-SQUID at fluxes 0, ±0.28$o and 
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FIG. 3; (a) The positions (/5o%) and (b) the widths (A/) 
of the histograms. Black solid and dotted lines are results 
(with known junction parameters) of TA and MQT model, 
respectively, ignoring dissipation. Blue, red and orange solid 
lines are the parameters of the simulated histograms based on 
the LO model discussed in the text. It assumes escape and 
possible relocalization events to be rare; thus it is not valid 
for data at $ = ±0.41$o, where the rate of phase difl^usion 
events approaches the relaxation in the well. 
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±0.41$o by horizontal dashed lines. It can be seen that 
the intersections of the dashed lines and the boundary of 
the UPD regime are very close to the experimental val- 
ues of Tq. The saturation of the histograms and their 
re-entrant steepness is thus a manifestation of the cross- 
over from TA escape into UPD due to dissipation. In the 
case of a single JJ with R{u!p) ~ 680 ft we obtain Q « 13 
yielding « 0.1 Ic and Tu « 1.2 K, instead of the mea- 
sured 0.3 Ic and 650 mK. By taking R{ujp) ~ 230 we 
obtain Q = 4.4, yielding /,„ « 0.3 Ic and Td « 700 mK 
instead, which is consistent with what we obtain from 
the position 1^1 1 c- 

Fig. 12 shows that the plain TA-MQT models cannot 
account for our observations. Except for the data at the 
lowest temperatures the width and the position of the 
measured histograms deviate from simulated ones (dot- 
ted line). Dissipation alone cannot explain the difference. 
The basic TA model yields A/ oc T^/^ ^ and it can 
be seen in Figs. El and El that the dc- SQUID has weaker 
temperature dependence even well below T^i. In the dc- 
SQUID there are just few energy levels in the well and 
thus the assumptions of continuous energy spectrum are 
not valid here ISj . The semiclassical model of Larkin and 
Ovchinnikov |3| takes into account the influence of the 
shape of the potential, in particular the fact that it is 
not harmonic. The total escape probability is calculated 
using Pesc{T) = 1 - /5fc(T), whcrc pkir) is the proba- 
bility of finding the particle in a state k after the current 
pulse of length r. The kinetic equation of the phase par- 
ticle can be written as ^ = Y^jilkjPj - IjkPk) - ^kPk- 
We take into account only transitions between neigh- 
boring levels and tunnelling out, Ffc. The relaxation 
rates between levels j and j ~ 1 are well approximated 
by 7j-ij- = ju!p/4:Q. We assume detailed balance and 
write 7jj_i=e~'^^^'~^^"'^^7j_i The positions and the 
escape rates are calculated using the results in Ref. 0. 
The final state p = [pi P2 ■ ■ ■ Pk] is calculated by nu- 
merically using p(r) — i e^'^*^p(0)(it, where A is the 
transition matrix. / is set to zero in the beginning, and 
p(0) is Boltzmann distributed. 

The effect of the rclocalizing dissipation can be taken 
into account also in the quantized energy level model. 
Using again Ed — 8Ej/Q, and the fact that the en- 
ergy difference between the two successive maxima is 
—2'kEjI/Ic: the level energy E must satisfy 

AU ~E < Ej (2nI/Ic - S/Q) (3) 

to allow switching into the free running state. If (jSjl is 
violated, the corresponding tunnelling rate is set to zero 
0. Note that Q gives the same threshold as JQ) for 
AU = E, but in Eq. ^ we assume that after tunnelling 
the starting point is not at the potential maximum. The 
solid red lines in Fig. Elpresent results of simulations with 
quantized energy levels and dissipation for the data of 
the dc-SQUID at zero and ±0.28*o fluxes. At ±0.41$o 
Ic is so small that the escape probability is large even 



at zero bias (except at the lowest temperatures). This 
means that the phase is moving constantly rather than 
infrequently escaping from a localized state, and thus our 
model does not work anymore. The only fitting parame- 
ter was Q, and the fitted values were in a very reasonable 
range. At <i> = we find Q « 6 at the lowest tempera- 
ture, and it decreases with increasing temperature up to 
4 at 325 mK. At ±0.28<i>o Q « 3 to 4, again decreasing 
with temperature. In Fig. Elis presented the 150% and AI 
parameters of the simulated histograms. We present also 
the simulated parameters for single junction {Q « 4) [loj . 
In the measurements on a single junction, the number of 
energy levels is large (« 20) and the level separation is 
smaller than the level width. Our model assumes well- 
separated levels and thus the agreement with a simple 
TA model is better for a single junction at temperatures 
T < Td, especially in Fig. El (a), whereas re-entrant 
steepness of the histograms could not be explained with 
the basic TA model. The agreement between the simula- 
tion and the measurements on a DC-SQUID is excellent. 
The position and the width of the measured histograms 
coincide and, in particular, at higher temperatures the 
simulated histograms for both samples also start to get 
steeper again. At higher temperatures the upper energy 
levels, whose escape rate is significant with smaller poten- 
tial tilting angles, are populated as well. The histograms 
thus peak at smaller currents and the condition (EJ is 
not necessarily fulfilled anymore. What remains in the 
measured (and simulated) histograms is the escape from 
the states at the tail of the Boltzmann distribution above 
the dissipation barrier. 
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